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Abstract. We study the average entanglement entropy of blocks of contiguous spins 
in aperiodic XXZ chains which possess an aperiodic singlet phase at least in a certain 
limit of the coupling ratios. In this phase, where the ground state constructed by a real 
space renormalization group method, consists (asymptotically) of independent singlet 
pairs, the average entanglement entropy is found to be a piecewise linear function 
of the block size. The enveloping curve of this function is growing logarithmically 
with the block size, with an effective central charge in front of the logarithm which is 
characteristic for the underlying aperiodic sequence. The aperiodic sequence producing 
the largest effective central charge is identified, and the latter is found to exceed 
the central charge of the corresponding homogeneous model. For marginal aperiodic 
modulations, numerical investigations performed for the XX model show a logarithmic 
dependence, as well, with an effective central charge varying continuously with the 
coupling ratio. 
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1. Introduction 

The entanglement properties of strongly correlated quantum systems is presently 
attracting great attention in condensed matter physics [H [21 El H] , quantum information 
theory [5], and DMRG theory [6]. From the point of view of condensed matter physics, 
one main aim is to understand the role of quantum correlations in quantum critical 
phenomena. A widely used estimator of quantum correlations is the entanglement 
entropy of a subsystem s embedded in a larger system, which is defined as 



where p s denotes the reduced density matrix of the subsystem. Near quantum critical 
points entanglement entropy was found to display interesting scaling behavior [U [2, H E] . 
In one dimensional homogeneous critical spin chains the entanglement entropy belonging 
to a block of L contiguous spins was shown to grow asymptotically as 



where c is the central charge of the associated conformal field theory, and A; is a 
non- universal constant, while Sl was found to remain bounded in non-critical chains. 
This type of scaling of the entropy seemed to be related to the conformal invariance 
of homogeneous critical systems [H [2] . However, similar logarithmic scaling of the 
entropy was also found, although, with a different "effective" central charge, in critical 
inhomogeneous systems such as systems with a single defect [21 EE] or quenched random 
disorder [31 HH [121 [131 EEU EES! ESI E] , where the ground state and the low temperature 
properties are quite different from those of homogeneous models |18j . Some of these 
recent results have already been used in the experimental identification of certain 
disordered phases [T9] . 

In the present paper a systematic study of the scaling of the mean entanglement 
entropy in a wide class of aperiodic models is carried out. The considered models, 
the self-similarly singlet-renormalizable aperiodic XXZ chains, are characterized by the 
property that the aperiodic sequence of the couplings remain invariant under a real 
space RG procedure which produces effective spin singlets. The ground states of these 
models are aperiodic singlet states, at least in the limit where the RG method applies. 
Some features of two particular models from this class (the Fibonacci and the tripling 
chain) have already been investigated in a recent work [2D] . 

We shall derive a general formula for the effective central charge characterizing 
the asymptotic L-dependence of the average entropy and show that in the XXZ chain 
this type of inhomogeneity may even enhance entanglement in contrast to single defect 
and random perturbations, in which cases the effective central charge decreases or 
remains invariant. Although, it is known that in some special chains of spins with 
many components also random disorder may increase entanglement [13] . 

The paper is organized as follows. In Secj2] some basic features of aperiodic 
sequences and the renormalization group method for inhomogeneous quantum chains 
are recapitulated and the class of models to be studied is introduced. In Secj3]we derive 




(1) 



S L =UnL + k, 



(2) 
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an analytical expression for the entanglement entropy asymptotics in the considered 
aperiodic singlet phases, whereas in SecJH results of numerical investigations for the XX 
chain are presented. Secj5] is devoted to the discussion of the results. Finally, in the 
Appendix, we prove that among the self-similarly singlet-renormalizable aperiodic XXZ 
models the one with strong tripling modulation has the largest effective central charge. 



2. Aperiodic quantum spin chains 

We study the antiferromagnetic spin-| XXZ chain defined by the Hamiltonian 

H = 22, Ji(SiSi+i + SfSl+i + A S*S* +1 ), (3) 

% 

where the Sf's [a = x,y,z) are spin-| operators, the site-dependent couplings J, > 
are taken from a finite set { J a , Jb, ■ ■ .} and are modulated according to some aperiodic 
sequence, whereas A (0 < A < 1) is the anisotropy parameter which is initially 
site- independent, but, as we will see later, it becomes site-dependent in the course of 
renormalization, as well. 



2.1. Aperiodic sequences 

The aperiodic sequences considered in this work are composed of letters taken from 
a finite alphabet {a, b, c, . . .} and are defined by an inflation rule which assigns to 
each letter a word (i.e. a finite sequence of letters) and the (infinitely many) repeated 
application of which generates the sequence. E.g. the well-known Fibonacci sequence 
consists of two different letters, a and b, and is defined by the inflation rule 

{a — > w a = ab , A s 

4 
o — > Wb = a. 

The first few inflation steps starting from letter a results in the following subsequent 
strings: a,ab,aba,abaab,abaababa, . . .. Many properties of aperiodic sequences are 
encoded in the substitution matrix M a/3 = n a (wp), where n a (wp) denotes the number 
of letters a in the word uip. E.g. the largest eigenvalue A + of M gives the asymptotic 
ratio of lengths of the subsequent strings in the inflation procedure. 

A finite aperiodic Hamiltonian is constructed by taking a finite section of the infinite 
aperiodic sequence and assigning a coupling J a to link + 1) whenever the ith letter 
of the sequence is a. 



2.2. Relevance/irrelevance criteria 

According to a heuristic criterion developed originally for random models [2TJ and later 
generalized to aperiodic systems [22]. the relevance of a weak aperiodic modulation is 
related to the wandering exponent u which characterizes the fluctuations of the local 
order parameter. To be concrete, an aperiodic perturbation is predicted to be relevant 
if 

»>l-l, (5) 
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where d is the dimension of the system and v is the correlation length exponent of the 
homogeneous system. 

For the particular case of aperiodic XX chains an exact relevance/irrelevance 
criterion is available which is obtained by an exact renormalization group method 
[23] . Here, the relevant wandering exponent is related to the fluctuations of the non- 
overlapping pairs of letters, and the aperiodic perturbation is relevant (irrelevant) if 
uj > (u < 0), in agreement with the Harris-Luck criterion with d — 1 and v — 1. For 
uj = 0, which is realized e.g. in the case of Fibonacci modulation, the perturbation is 
strictly marginal and critical exponents vary continuously with the coupling ratio. 

According to a weak perturbation renormalization group study of the general XXZ 
chain [23] the weak Fibonacci modulation drives the system away from the pure system 
fixed point for any anisotropy < Aq < 1. 



2.3. Strong disorder renormalization group 

Recently, many results have been obtained for aperiodic quantum spin chains by means 
of a real space renormalization group method J25J [26] which was originally developed 
for random systems [271128]. The main point of the method is the successive elimination 
of high-energy degrees of freedom connected to the largest coupling in the system. 
For the antiferromagnetic XXZ chain, the spin pair (Si,Si+i) coupled by the largest 
bond Jj, which is approximatively in the singlet state ^=(| f j) — | if)) provided that 
Ji Ji-i, Ji+i, AjJj ^> Aj_i Jj_i, Aj + i Jj+i, is eliminated from the chain, and effective 
new couplings between spins SV-i an d S i+2 

~ Jj_iJj + i ~ 1 + Aj 

J = (TTA^V A = -2-^-1^ (6) 
are generated by second order perturbation theory (see Fig. [I]). This elimination step 

.... — ► -h... r\ ..^ 

Jj-l ^ Ji+1 

J 

Figure 1. Renormalization scheme in the XXZ chain. The spins connected by the 
strongest bond Ji form a singlet and a new bond J is generated. 



is iterated, leading to the gradual decrease of the number of active spins and that of the 
energy scale set by the actually largest effective coupling. For the critical random XXZ 
chain the distribution of couplings is broadening without limits during the procedure, 
implying that the conditions for the perturbation calculation are better and better 
fulfilled, and it becomes asymptotically exact in the (infinite randomness) fixed point. 
The latter is attractive for any amount of disorder, in agreement with the Harris criterion 
which predicts the random inhomogeneity to be relevant for < Ao < 1. The resulting 
ground state constructed by the method consists of singlets of spins which are separated 
by arbitrarily large distances, and the system is in the so called random singlet phase. 
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The method was applied to the XXZ chain with various types of aperiodic 
modulations, where identical blocks are simultaneously renormalized in the whole chain. 
After such a renormalization step the structure of certain aperiodic modulations remain 
self-similar (25| [26] . E.g. for the Fibonacci XXZ chain which is known to be such 
a system, the blocks to be renormalized correspond to the words aba and ababa (see 
Fig. [2]), and the renormalization step is essentially the reversed process of the triple 
application of o F given in eq. (jl]) |25|. When A = 1, the ratio of the couplings 

• • • • • • • • • 

Ja Jb Ja Ja Jb Ja Jb Ja 

T T 

... r\ ... r\ r\ ... 

Jb Ja 

Figure 2. Renormalization of the Fibonacci XXZ chain. 

r a b = J a / Jb and the anisotropy renormalize to r* ah = and A* = 1, respectively. 
The method becomes asymptotically exact in the fixed point which is attractive for 
any r ao < 1 in agreement with the fact that weak perturbation is relevant. The 
ground state consists of spin singlets the possible lengths of which form a subset of the 
Fibonacci numbers and the system is in the so called aperiodic singlet phase. When 
A = 0, the ratio r a b remains invariant in the RG procedure, indicating that the 
Fibonacci modulation is marginal. Here, the decimation steps are exact only in the 
strong modulation limit — ► 0. For the case < A < 1 the coupling ratio and 
anisotropy parameter renormalize to some r* ab > and A* = 0, respectively. Thus, 
the strong perturbation is expected to be marginal whereas the weak perturbation is 
relevant. Here, the critical exponents are non- universal again, depending both on r ao 
and A [25J. 

2.4- Singlet producing sequences 

In this work, we focus on XXZ chains with such aperiodic sequences (like the Fibonacci 
chain), where the RG procedure produces exclusively effective singlets and the sequence 
of couplings remains invariant. Generally, in such systems, the blocks to be renormalized 
consist of an even number of spins (an odd number of couplings), and the couplings in 
each block are alternately weak and strong (with strong bonds at even places) in order 
to avoid the forming of effective spins [26J. Furthermore, not only the RG-invariance of 
the sequence of couplings (self-similarity) is required, but also the order of strength of 
the couplings should remain invariant under renormalization, i.e. J a > Jp should imply 
J a > Jp- Each such sequence has a unique canonical inflation rule where the word 
w a represents that block which just renormalizes to the two-spin block represented by 
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letter a, that is, the inflation rule just corresponds to the reversed renormalization 
step. To avoid some unphysical situations we also require two technical conditions: the 
ergodicity of the inflation rule, i.e. the same infinite sequence should be obtained in 
the bulk independently of the letter which the iteration is started from, and the RG 
procedure should be "unambiguous", i.e. w a ^ w@ if a ^ [3. We call the sequences 
fulfilling the above conditions singlet producing self-similar sequences. 

The family of singlet producing self-similar two-letter sequences have the following 
canonical inflation rule: 

{a — ► w a = a6a(6a) m_1 , . 

b - w b = a(ba) n -\ [ ) 

where the integers m, n satisfy 1 < n < m, and letter b now and in the rest of this work 
represents the strongest bond J - The sequence generated by the transformation with 
m = n = 1 is the well-known silver-mean sequence, the choice m = n = 2 generates the 
Fibonacci sequence, whereas the transformation with m = n = 3 can be expressed as 
a composition of simpler inflation transformations as o~p ° &n ° <Jf- An example for a 
singlet producing sequence composed of three letters is 

!a — > w a = abababc 
b — > Wb = abc (8) 
c — > w c = ababc, 

which is the sequence obtained after the first RG step of the 6-3 sequence [26]. Another 
one is the tripling sequence 

!a — ► w a = aba 
b — ► Wb = cbc (9) 
c — > w c = abc. 

which plays a special role among singlet producing sequences as will be shown later. 
For the latter two sequences the relations Jb > J c > J a remain invariant under the 
renormalization. Generally, the new effective couplings after an RG step are given in 
terms of the couplings of the previous generation and the canonical substitution matrix 

as 



3 rx 



Y[A^ Wa \ (10) 

which are obtained by eliminating the strong bonds Jb one by one in the blocks 
represented by the words w a , using eq. ©• The considered sequences can be categorized 
into three types on the basis of their relevance in the strong perturbation limit. 

Those sequences for which the strong perturbation is relevant independently of 
Ao are of type I. In this case, the coupling ratios renormalize to zero, the procedure 
is asymptotically exact and the resulting ground state is an aperiodic (independent) 
singlet state. These aperiodic modulations are exactly known to be relevant in the XX 
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model for arbitrary coupling ratios, and this is expected to hold also for < A < 1. 
The two-letter sequences with m > n and e.g. the 6-3 sequence are of type I. Type 
II sequences, such as the Fibonacci sequence, are relevant for the Heisenberg chain 
(A = 1), but are strictly marginal for the XX chain (A = 0). Here, for A = 0, the 
coupling ratios remain invariant under the renormalization, whereas for intermediate 
initial anisotropics < A < 1 the coupling ratios tend to a finite limiting value. So, 
the ground state is an aperiodic singlet state for A = 1, but for A < 1 this is true only 
in the limit of vanishing initial coupling ratios. The two-letter sequences with m = n 
are of this type. Type III modulations are strictly marginal for the XX model, and 
the coupling ratios remain invariant under renormalization, independently of Ao, thus 
marginality is indicated (at least in the strong modulation limit). Here, the ground state 
is an aperiodic singlet state only in the limit of vanishing initial coupling ratios. The 
tripling sequence and other sequences with words of equal length belong to this class. 

3. Entanglement entropy in the aperiodic singlet phase 

3.1. Analytical expression for the scaling of entanglement entropy 

In this section, we consider an infinite chain and calculate the average entanglement 
entropy of subsystems consisting of L consecutive spins. In the aperiodic singlet phase 
this can be done analytically, since, similarly to the random singlet phase, one simply 
has to count the singlet bonds connecting the subsystem with the rest of the system 
[3]. The entanglement entropy of one spin in the singlet pair is In 2, therefore each such 
bond gives a contribution of In 2 to the entanglement entropy (Fig. [3]). At any stage of 



the renormalization procedure, the spin pairs, which form singlets, are those connected 
by the actually strongest (and, at the same time, shortest) bonds, Jb. The evolution 
of bond lengths during the renormalization can be determined by arranging them in a 
vector w and noticing that the corresponding bond lengths after a renormalization step 
are obtained by multiplying w by the transpose of the canonical substitution matrix M. 
The length of the singlets in the nth generation (i.e. those generated in the nth step of 
the renormalization), which is equal to the length of the bond corresponding to letter b 
before the nth step, can thus be calculated after diagonalizing M, yielding 




L 



Figure 3. Calculation of the entanglement entropy. In the present case, the 
entanglement entropy of the block indicated by the arrow is 3 In 2. 
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where A + is the largest eigenvalue of M and u a denotes the component a of the left 
eigenvector of M associated with the eigenvalue A+. This eigenvector and the right 
eigenvector v of M associated with A + are normalized as follows 

J2 v » = J2 u " v " = l - ( 12 ) 

a a 

The terms in the r.h.s. of eq. (fill) , denoted by 0(A™ _1 ), involve the subleading 
eigenvalues of M and are negligible compared to the leading order term for large n. 
The concentration of singlets of the nth generation, i.e. their average number per unit 
length, is given by 

Pn = V b X+ {n ~ 1] . (13) 

The fraction of the chain which is covered by singlets of the nth generation is w n = l n p n , 
which tends rapidly to u b v b for large n. This quantity gives the probability for breaking 
a singlet bond of the nth generation when the chain is cut at a randomly chosen link. 
Therefore, each generation for which l n < L contributes a term of 2w n In 2 to the average 
entanglement entropy. The total contribution from singlets with l n < L is thus for large 
L 

Si(L) « 2n(L)u b v b \n2 + const, (14) 
where n(L) is the number of generations with l n < L. From eq. (fTTj) we get for large L 

ln(L/w 6 ) 



n(L) 



ln A. 



(15) 

- int 

where [x]i n t stands for the integer part of x. The contribution of singlets with l n > L is 

oo 

S 2 {L)= 2p n Lln2 = LA; n(L)_1 ln2, (16) 

n=n(L)+l 

where we have used the relation 

OO 1 

£Pn = -, (17) 

71=1 Z 

which expresses the fact that all spins are paired in the aperiodic singlet phase. Finally, 
we obtain for the asymptotic L-dependence of the entanglement entropy: 

S{L) = 2n(L)u b v b ln 2 + LA;" (LM ln 2 + const. (18) 

As opposed to translation invariant systems, S(L) is a piecewise linear function with 
breaking points at all possible singlet lengths l n , see Fig. HI The appearance of the 
characteristic lengths l n in the correlation functions is common in the aperiodic singlet 
phase Nevertheless, the enveloping curve of S(L), which asymptotically fits to the 
breaking points, grows logarithmically with L, hence one can write the entropy in the 
form 

S{L) = C -f\nL + P{L), (19) 
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where P(L) is an oscillating log-periodic function. From eqs. (j!5p and (ITS]) we obtain 
the general expression for the effective central charge 

c cff = 61n2-^L, (20) 
in A + 

which differs from the central charge of the homogeneous system c = 1 and it is 

characteristic for the underlying aperiodic sequence. Evaluating eq. ( 1201) for the two- 




300 



Figure 4. Average entanglement entropy as a function of the block size L given in 
cq. I|18p. evaluated for the Fibonacci sequence (solid curve). The dashed curve is the 
enveloping function S env (L) = £rf£ InL + const. 



letter sequences a mn yields 

61n2 (A + -m-l) 2 

Ceff - ; 7 —77 TTo, (/I J 

m A + mn + (A + — m — l) z 

with the rescaling factor A + = |(m + n + \J (m + n) 2 + 4(m — + 1)). For the 6-3 
sequence we get c c s = flrf ~ 0.6201, whereas for the tripling sequence c e fr = x^-f ~ 
1.2618. 



m=n=l (silver mean) 


0.6910... 


m=n=2 (Fibonacci) 


0.7962... 


m=n=3 


0.7819... 


m=n=4 


0.7519... 


m=n=5 


0.7228... 



Table 1. Effective central charges for aperiodic sequences defined in eq. ((7]) in the 
aperiodic singlet phase. 
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abaababaababaabaababaababaaba 





/rf^ ^ ^ ^ ^ /rf^i 

abaabaaabaabaaabaabaaabaaba 





abacbcabcabccbcabcabacbcabc 



Figure 5. The singlet bond topology of the Fibonacci (upper figure), the silver- mean 
(middle figure), and the tripling (lower figure) sequences. 



3.2. The maximally entangled aperiodic singlet topology 

Next, we turn to discuss the possible bounds of the effective central charge. Considering 
the two-letter sequences, a mn , for the sake of simplicity with m = n, we have found that 
after the maximum at the Fibonacci sequence (m = n = 2) the effective central charge 
is decreasing with the length of the words (see Table I.) and it finally tends to zero 
in the limit n = m —>■ oo. This is easy to understand, since the words themselves 
are periodic, and the blocks corresponding to the words have a dimerized (non-critical) 
ground state. Thus, the effective central charge can be arbitrarily close to zero. One 
may ask whether sequences (belonging to the considered class) with an arbitrarily large 
effective central charge can be constructed. The answer is negative. To see this at a 
heuristic level, note that the expression for the effective central charge in ( 1201) may be 
separated into two (non-independent) factors, the covered fraction and the rescaling 
factor. When the singlets of the same generation lie densely, leaving only a small 
concentration of unpaired spins between them, the covered fraction will be large which 
is favorable for strong entanglement, however, the rescaling factor will also be large, 
which has the opposite effect. Conversely, for a large concentration of lonely spins, 
the rescaling factor is moderate, however, the covered fraction is also small. Thus, one 
expects that there may be an optimal aperiodic singlet topology which maximizes the 
entanglement entropy, and indeed this is the case. We can prove that the aperiodic state 
belonging to the tripling sequence has the highest central charge (see the Appendix). 
Beside this proof, one can also give a less formal, but more suggestive argument for this. 
The singlet arrangement produced by the tripling sequence (see Fig. [5]) is optimal in 
the following sense: Regarding those subsequent renormalization steps which generate 
non-overlapping singlets (e.g. for the silver mean sequence pairs of subsequent steps), as 
one single renormalization step (in which singlets of different lengths may be produced), 
the number of unpaired spins between singlets of the same generation is at most one. 
For the tripling sequence there is an unpaired spin between every neighboring singlet, 
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4. Numerical analysis 

In order to check the results obtained in the previous section and to get some insight into 
the behavior of entanglement entropy for marginal perturbations when the ground state 
is not an independent singlet state and the RG breaks down, we performed numerical 
calculations for the XX model (A = 0). This model can be mapped to a system of free 
fermions [29] and the reduced density matrix of the subsystem can be written as the 
exponential of a free-fermion operator [30]. The entanglement entropy can be computed 
from the eigenvalues Aj of the correlation matrix of fermion operators CV,- = (c\cj) 
restricted to the subsystem by the formula [TJ [30] 

S(L) = -£[A* In Aj + (1 - A,) ln(l - A,)] . (22) 

i 

The computation of the correlation matrix necessitates the diagonalization of an N x N 
matrix, then diagonalizing the restricted correlation matrix of size L x L one obtains 
the von Neumann entropy from eq. ( |22l) . 

We performed numerical calculations for chains with both free and periodic 
boundary conditions, and calculated the von Neumann entropy of blocks which were 
located in the former case in the middle of the chain. The size of blocks varied from 2 
to N/2 and the largest system size we considered was N = 1024. The entropies were 
then averaged over a few thousand independent samples, whereas for periodic chains an 
averaging over the different starting positions of the block was also carried out. 



1.8 
1.6 

5 1 - 4 

CO 

1.2 
1 

0.8 

1 2 3 4 5 6 7 
InL 

Figure 6. Numerically calculated average entanglement entropy as a function of the 
block size L for the relevant two- letter sequence defined in (TTJ) with m — 2 and n = 1. 
The calculations were performed with free boundary condition. The size of the system 
is N = 1024 and the data were averaged over 10 4 samples. The solid line is the 
analytical prediction (fl~8|) . For this sequence c c g « 0.4459. 
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1 2 3 4 5 6 
InL 



Figure 7. Numerically calculated average entanglement entropy as a function of the 
block size L in the Fibonacci XX chain, for different coupling ratios. The calculation 
was performed with periodic boundary condition. The size of the system is N = 1024 
and the data were averaged over 20 independent samples and over the different starting 
positions of the blocks. The solid line is the asymptotics of S(L) for the homogeneous 
case r a b — 1, whereas the dashed curve is the expected asymptotical curve in the limit 
fab — > 0, given in eq. (HI). 

The simplest sequence among the two-letter sequences a mn which is of type I and is 
thus relevant even in the XX limit is the sequence with m = 2 and n = 1. The average 
entropy as a function of the block size, shown in Fig. El tends to the asymptotic curve 
characteristic for the independent singlet phase after a crossover region even for a finite 
initial coupling ratio in accordance with the expectations. 

For type II sequences, which are marginal in the XX limit and relevant in the 
Heisenberg limit, one expects that, analogously to the critical exponents, the effective 
central charge varies with the coupling ratio and the anisotropy parameter for Ao < 1. 
Numerical results in the XX case for the size dependence of the entanglement entropy for 
one representant of this class, the Fibonacci sequence, are shown in Fig. [7] for different 
coupling ratios r ab . As expected, the effective central charge, which is related to the slope 
of the straight line fitted to the breaking points of the curves in the semi-logarithmic 
plot, varies continuously with the coupling ratio, see also the numerical results in [20J. 
For an other type II sequence, the silver-mean sequence, the effective central charge as 
a function of the coupling ratio is shown in Fig. [HJ 

Type III sequences, being marginal for any A , are also expected to have effective 
central charges which may depend on the coupling ratios and the anisotropy. Numerical 
results for the tripling sequence in the XX limit, presented in Fig. [91 are in agreement 
with these expectations. 
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r ab 

Figure 8. Effective central charge for different coupling ratios r a b calculated for the 
silver-mean XX chain with periodic boundary condition. The arrow indicates the 
exactly calculated value in the aperiodic singlet phase, taken from TableHJ Inset: The 
average entropy for different block sizes L = N/2 which are chosen to be the first few 
singlet bond lengths l n in order to avoid log-periodic oscillations. The averaging was 
performed on the L different starting positions of the block. 




1 2 3 4 5 6 7 
InL 



Figure 9. Numerically calculated entanglement entropy as a function of the block size 
L in the tripling XX chain, for different values of the coupling ratios J a / J c = J c / Jb = r. 
The calculations were performed with free boundary condition. The size of the system 
is N — 1024 and the data were averaged over 10 4 samples. The solid line is the 
asymptotics of S(L) for the homogeneous case r = 1, whereas the dashed curve is the 
expected asymptotical curve in the limit r — ► 0, given in (fT8l) . 

5. Discussion 

We have studied in this work the average entanglement entropy of blocks of contiguous 
spins in a class of aperiodic quantum spin chains, where the sequence of couplings 
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remains invariant under the strong disorder RG procedure which produces effective 
spin singlets. For these systems, the ground state is, in the limit where the RG method 
applies, an aperiodic (independent) singlet state. Knowing the topology of singlet bonds, 
the average von Neumann entropy as a function of the block size can be calculated 
analytically. This quantity, apart from the log-periodic oscillations, is found to grow 
logarithmically with the size of the block. A general expression for the effective central 
charge was derived, and the modulation producing the largest central charge was found. 

Beside the singlet producing self-similar sequences studied in this work, which 
can be termed as "pure", also other singlet producing sequences can be constructed 
by composing canonical inflation transformations of different pure sequences. An 
example is the transformation o~ n o a 33 which is, on the other hand, just the 
square of the transformation cfs = an o o~ F) therefore both generate the same 
"golden-silver- mean" sequence. For composite sequences involving m consecutive pure 
transformations, M 1; M 2 , . . . , M m , the calculations outlined in Sec. [3] can be generalized 
in a straightforward way. In this case the resulting effective central charge is the 
sum of central charges calculated from the m cyclic permutations of the product 
Mi o M2 o . . . o M m using eq. (1201) . For example, for the golden-silver-mean sequence we 
obtain c cff = ^ + - 4 % ln2^ 0.6673. 

Finally, we mention that for sequences of type I and II, which are relevant at least 
for Ao = 1, where the perturbation of arbitrary strength is expected to drive the system 
to the strongly aperiodic fixed point, we have not found larger values of c e fj than that of 
the homogeneous model c = 1. This result is consistent with the "generalized c-theorem" 
proposed in Ref [3j as a possible extension of the original c-theorem [31] to disordered 
systems. According to this conjecture, the effective central charge, defined through 
entanglement scaling, decreases both for homogeneous and disordered systems along 
RG trajectories induced by relevant perturbations. The sequences for which we have 
found c e fi > 1, the tripling and a sequence with words of length five, both belong to type 
III. These types of perturbations are marginal at least at the strongly aperiodic fixed 
point (in the special case Ao = even for perturbations of arbitrary strength), therefore 
these cases are out of the range of validity of the generalized c-theorem. However, for 
Ao > 0, the behavior of the effective central charge for finite perturbations of this type 
is still an open question. It is worth mentioning that our results show that disordered 
marginal perturbations can increase the effective central charge, which is not expected 
in case of homogeneous marginal perturbations [32] . 
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Appendix A. 



Here, we prove that the tripling state has the largest effective central charge among the 
states generated by singlet producing self-similar sequences. 

First, we notice that the relations ( fl3l) and (fT7j) imply i> 6 = |(1 — A^ 1 ). Inserting 
this into the expression of the central charge in eq. f[2"Ul) we obtain 

(1 " Ajl) "i. (A.i) 



c cS = 3 In 2- 



lnA 



+ 



In the considered class of models, as we have mentioned, the strongest effective bond is at 
the same time the shortest one (or, more precisely, its length is smaller than or equal to 
that of the other effective bonds) after each renormalization step. Therefore, the length 
l n+ i of singlets produced in the n + 1st renormalization step cannot exceed the average 
distance £ n between active (unpaired) spins after the nth step, i.e. £ n > l n+ \. The 
concentration of active spins decreases exactly by a factor of A + in each step, therefore 
£ n = A n , which, together with eq. ffTTl) leads to w b < 1. Putting this in eq. (1A.1I) yields 

(A.2) 



c e ff < 31n2- 



lnA_ 



Note that we have immediately an upper bound for all the considered sequences, since 
the r.h.s. of (1A.2j) is monotonously decreasing with A + (> 1), hence c e g < 3 In 2. 
However, this bound can be further sharpened, as we show below. 

For A + > 3, the statement we want to prove follows directly, since in this case the 
maximal value of the upper bound ( 1A.2I) is |^|, which is just the effective central charge 
of the tripling state. 

If A + < 3, then letter b which represents the strongest bond must be transformed by 
the canonical inflation rule into a single-letter word. Indeed, Wb = abc . . . would imply 
l b n+l = l* + l b n + l c n + . . . > 3Z£, which is equivalent to A+ > 3 . Here, l a n = u a \\ + C(A™ ) 
denotes the length of the bond represented by letter a (or shortly, bond a) after the nth 
renormalization step. The most general form of canonical inflation rules with A + < 3, 
is thus 



a {1) 
a m-l 

«£ } 

aS 2) 
a?) 



(2) 
,.(2) 



"2 



J 1 ) 



(2) 
(2) 



w (2) 
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(k) 
a\ 

(k) 



(k) 

a n k ~l 

a™ 



(fc) 



w 



,(*) 

(*) 



where denotes a single letter, specially ai = b is the letter representing the strongest 
bond, whereas stands for a word comprising at least three letters. The indices n.j 
are arbitrary positive integers with the restriction % > 1. In the following, we use the 
notation for the density of letter and write A™ for the asymptotical length of 
the corresponding bond after the nth renormalization step. With these notations, eq. 
flTJ) reads 



i=l 4=1 

Furthermore, we need the following relations 



(A.3) 



J i+1 



u 



0) _ \ . (J) 

i+l — A + U i 5 



(A.4) 
(A.5) 

A+ > V3, (A.6) 

the first two of which follow simply from the form of the inflation rules. Applying eq. 
(1A.5j) rii — 1 times for j = 1 yields 

<W - A!TV (1) . (A.7) 



ni 



"1 • 



Using that a$ is transformed to the long word and that v>i < uj 3 '' for all z, j, we 
have A + u^ > 3W] 1 . This inequality together with eq. (I A. 71) give immediately relation 
(TAJ]). 

Now, we consider the group of letters with a fixed index (j) and introduce 



the quantities v^> = 
relations flA.5j) imply 



Yh=i W and u> 



i j) /v ij \ such that E^i 



1 < 



re 



U) 



w{ j) 



(?) 



w 



U) 



Ai< 



< 



n, ,n 3 — 1 



(?)' 



1 holds and 



(A.8) 



With these variables we can write 



£4 



S J) £ ^ A 1 " 1 = u^t;^ £ wfA 



(A.9) 



i=l i=l i=l 

where the first equality follows from eq. (1A.5[) . The latter expression can be written 
also in the form 
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One can show that, since A + > 1, a lower bound of this expression is obtained by 
substituting the relations flA.8j) as equalities into the fraction on the r.h.s. of eq. flA.lOj) . 
and we get 

5>{V> > u^\ ^ X \^ X+ - l) . (A.11) 

i=l A + - 1 

We know that > u[ , however, we can sharpen this bound if rij < m, where m is the 
positive integer for which a/3 < A + < m ^/3. (From (1A.6j) it follows that m > n\ > 1.) 
In this case, ( 1A.5[) yields > /A+. Using these relations we obtain 

Tlj 

E^ ) ^ ) >nS 1 V^P,(A + ), (A.12) 
i=i 

with 



" J ,' (A+ " 1) if n,- > m 

J to — - if n,- < m. 

a+(a + j -i) J 

Summing over the index (j) and using the normalization ( 1A.3I) we get 

1 > mS X) E vU)p j( x +) > u( i ] min {^(A+)}, (A.14) 
i=i 

where the latter inequality follows from Y^j=x^ — !• It is easy to check that in the 
interval a/3 < A + < m ~\/3 

(1-AT 1 ) 2 , . 

< — — (A. 15) 



P i (A+)lnA+ ~ 3 In 3 

holds for all Pj(X + ). From this inequality and (1A.14I) we obtain an upper bound on Ui 

— f^sid^-)- (A - 16) 

Substituting this into the formula of the effective central charge (1A.1I) yields 

In 4 , , 

«- £ hll- (A17) 
The obtained upper bound is again just the value of the effective central charge of the 
tripling state. Thus, we have completed the proof. 
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